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Abstract:  Cylindrical  plasma  pinch  experiments 
are  conveniently  done  in  coaxial  geometry.  However, 
when  a  coaxial  line  is  driven  by  a  nonaxi symmetric 
source  of  current,  there  result  asymmetries  in  the 
magnetic  field  in  the  coaxial  line  which  in  turn 
will  drive  Rayleigh-Taylor  instabilities.  In  this 
paper  we  calculate  the  field  asymmetry. 

Introduction 

Our  basic  goal  is  to  vaporise  and  then  implode 
a  thin  cylindrical  metal  foil  by  means  of  a  very 
large,  rapidly  rising,  current  pulse.  The  resulting 
shell  of  plasma  is  compressed  as  a  z-pinch.  Because 
of  the  geometry  of  the  foil  it  is  natural  to  feed  it 
by  means  of  a  coaxial  line.  The  foil  then  simply 
forms  a  section  of  the  inner  conductor  of  this  line. 

In  ideal  coaxial  geometry,  the  z-current 
(current  parallel  to  the  axis  of  the  coaxial  line) 
and  the  resulting  Bg-field  are  both  perfectly 
axially  symmetric.  At  present,  however,  our  current 
sources  are  parallel  plate  generators.  Connecting 
such  a  generator  to  a  coaxial  line  leads  of 
necessity  to  a  departure  from  axial  symmetry  of  the 
whole  system.  The  result  is  that  the  current 
density  and  the  magnetic  field  at  the  foil  are  no 
longer  axisymmetric  and  this  will  drive  a 
Rayleigh-Taylor  instability  (a  flute  mode)  of  the 
plasma  shell.  Such  an  instability  could,  if  large 
enough,  severely  limit  the  maximum  compression  and 
temperature  achieved  in  the  pinch.  Before  the 
Rayleigh-Taylor  problem  can  be  undertaken,  we  must 
know  the  details  of  the  driving  pressure  asymmetry, 
and  this  is  the  subject  of  this  paper. 

The  Magnetic  Field  in  the  Dielectric 


$  =  Z  [Am(r/r0)Dt  -  Bra(r0/r)m]  sin  m6  ,  (1.3) 
m=l 

where  A  and  B  are  arbitrary  functions  of  time 
mm 

alone  and  r^  is  a  typical  (constant)  radius, 

introduced  for  convenience.  Using  Eq.  (1.3),  the 
magnetic  field  components  are 

Br  =  1  7  lAm(r/r0)m  +  Bm(r0/r)"]  Si"  m6  ’ 

m=l 

B0  =  2I/r 

+  Z,  7  tAm(r/r0)nl  '  Wr)>1  C0S  K 

m=l  (1.4) 

B  =  0 
z 

between  the  parallel  plates. 

Between  a  pair  of  coaxial  cylinders  r  =  a  and 
r  =  b  (b  =  a+h) ,  and  assuming  the  same  angular 
symmetry  as  in  Eq.  (1.3),  we  may  write  the  solution 
as 

oo 

*  =  Z  [cV^2  -  D“e-A-Z]z  (Ar)  sin  m6  ,  (1.5) 

,  m  m  m 

m=l 

where  Z  (Ar)  is  an  arbitrary  linear  combination  of 

m  iV  * 

the  Bessel  functions  J  and  Y  and  C  ,  D  are 

m  m  m  m 

arbitrary  functions  of  time.  Strictly  speaking 
Eq.  (1.5)  should  also  include  a  summation  over  A.;  we 
shall  take  care  of  that  presently.  The  boundary 

condition  that  B  vanish  at  r  =  a,  the  inner 
r 

cylinder  wall,  is  taken  care  of  by  choosing 


In  our  problem  it  will  be  convenient  both  in 
the  parallel  plate  sections  and  in  the  coaxial 
sections,  to  work  in  cylinder  coordinates  (r,  0,  z). 
Then  the  plate  conductors  will  have  the  equation 
z  =  a  constant  whereas  on  cylindrical  conductors 
r  =  a  constant.  We  shall  represent  the  magnetic 
field  in  the  dielectric  region  between  conductors  as 


3  =  2I/r  efl  +  V0  ,  (1.1) 
where  e^  is  a  unit  vector  in  the  0-direction  and  I 

is  the  total  z-current  flowing  inside  the  radius  r, 
and  <)>  is  a  potential  function  describing  the 
asymmetry.  We  have  neglected  the  displacement 
current,  i.e.,  the  capacitative  impedance.  We  shall 
not,  however,  assume  that  9B/9t  =0.  As  the 
divergence  of  B  must  vanish,  we  have 


rar^P  r2ae2  9z2 


(1.2) 


This,  plus  appropriate  boundary  conditions,  will 
determine  $  and  thence  B. 

First  let  us  consider  the  region  between  two 
parallel  plates,  the  upper  located  at  z  =  z^  and  the 
lower  at  z  =  zQ  -  h.  The  dielectric  thickness  h 
will  later  be  assumed  to  be  small.  The  boundary 
condition  on  the  magnetic  field,  namely  that  B  =  0 


at  z  =  Zq  and  at 
by  the  solution 


h,  is  trivially  satisfied 


Z  (Ar)  =  Y'(\a)  J  (Ar)  -  J' (Aa)  Y(Ar)  ;  (1.6) 
m  m  m  m  m 


then  Z'(Aa)  =  0.  Here  primes  mean  derivatives  with 
m 

respect  to  the  argument,  Ar. 

Next  we  must  satisfy  the  boundary  condition 

that  B  vanish  at  r  =  b  =  a+h.  It  is  here  that  we 
r 

make  use  of  the  fact  that  h  is  small  by  Taylor 
expanding, 

Z' (Ab)  =  Z' (Aa)  +  Ah  Z"(Aa) 
mm  m 


+  ^(Ah)2  Z" ' (Aa)  +  •••  .  (1.7) 

m 

Now  we  can  evaluate  the  derivates  on  the  right  hand 
side  by  noting  that  Zm(p)  satisfies  Bessel's  equa¬ 
tion  and  that  Z' (Aa)  vanishes  by  its  construction, 
m 

In  this  fashion  we  obtain  the  expansion 

(  2 


Z;(Ab)  =  ZjAa) 


(-f~2  -  D(*h) 

\  a 


+ 


1 

2Aa 


(-f^  -  1) (Ah)2  +  •••( 
A  a  ) 


(1.8) 


When  Ah  is  sufficiently  small,  we  can  break  off  the 
sum  in  the  curly  brackets  after  a  modest  number  of 
terms  and  satisfy  the  boundary  conditions  Z' (Ab)  =  0 
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by  equating  this  bracketted  expression  to  zero.  To 
the  order  actually  written  out,  this  determines  A  as 

A  =  (m/a)  [1  -  %(h/a)  +  0(h2/a2)]  .  (1.9) 

What  is  required  for  this  to  work  is  that  (h/a)  must 
be  small.  If  we  define  the  mean  radius  p  as 

p  =  a(l  +  |^)  =  h(a  +  b)  ,  (1.10) 

then  Eq.  (1.9)  can  be  written 

A  =  m/p  (1.11) 

and  Eq.  (1.8)  states  that  Z^(mb/p)  =  0(m[h/p]  ),  so 

our  boundary  condition  is  satisfied  and  for  all 
b=a+h>r>a,  Z  (mr/p)  =  Z  (m) ,  both  to  high 
precision.  We  can,  therefore,  replace  the  radial 
functions  by  constants  and  define  new  arbitrary 
functions  which  depend  only  on  time 

C  (t)  =  C"Z  (m)  ,  D  (t)  =  d“z  (m)  .  (1.12) 

m  mm’  m  mm 

The  fields  are  then  given  by 


B  =  21/ r  +  2  (m/r)  [C  emz/P 

o  .  m 

m=l 

-mz/p 

-  D  e  ]  cos  m0  ,  (1.13) 

m 

B  =  Z  (m/p)  [C  emz/P  +  D  e  mz/P]  sin  m0 
z  ,  r  1  m  m 

m=l 

to  high  precision. 

Boundary  Conditions  at  the  Junction  of  Two  Lines 


In  Fig.  1  we  have  sketched  the  junction  of  a 
parallel  plate  line  with  a  coaxial  line. 


To  work  out  the  details  of  the  transition 
region  we  would  need  to  find  expressions  like 
Eqs.  (1.4)  and  (1.13)  for  this  region  and  match  the 
fields  smoothly.  We  shall  not  do  this  but  shall 
rather  make  use  of  the  fact  that  this  troublesome 
region  is  small  (h/a  «  1). 

First  of  all,  we  can  equate  the  two  expressions 

for  B„  evaluated  at  the  corner  z  =  zA,  r  =  b 
o  U 


obtaining 


m  mz./p 

B  (r./b)m  =  C  e  -  De 

m  0  m  m 


-mz0/p 


(2.1) 

We  need  one  more  boundary  condition  which  we 

shall  express  as  flux  conservation.  First  of  all  we 

consider  the  surfaces  Z  and  2  to  be  extended  a 

p  c 

distance  6,  the  skin  depth,  into  the  metal  walls. 
The  surface  areas  included  within  the  small  angle 
d8,  after  this  extension,  will  be 


dZ  =  b(h  +  26)  d0  , 

P 

dZ  =  3g[(b  +  6)2  -  (a  -  S)2]  d6 
c 


(2.2) 


where  p  is  the  mean  radius  given  by  Eq.  (1.10).  We 
have  assumed  that  the  skin  thickness  is  small  enough 
that  it  is  not  influenced  by  the  curvature  of  the 
cylindrical  surfaces.  Our  flux  conservation  condi¬ 
tion  becomes 


-B(p£)  dZ  =  B(cyi:)  d2  .  (2.3) 
r  p  z  c 

The  minus  sign  arrises  because  an  inward  (i.e., 
negative)  radial  flux  continues  as  an  upward  (i.e., 
positive)  z-flux.  From  Eqs.  (1.4),  (1.13)  and  (2.2) 
we  see  that  Eq.  (2.3)  becomes 


mz„/p  -mz_/p 

+  Bm('r0/b)in  =  _Cme  ”  Dme 

m  U  m  U  m  m 

(2.4) 


Note  that  the  finite  skin  thickness  cancelled  out 
exactly. 

The  boundary  conditions  (2.1)  and  (2.4)  give 


A 

m 


=  -Dm(rn/b) 

m  0 


m 


-mz0/p 


B 

m 


-C  (b/r  )m 
m  0 


mz0/p 

e 


(2.5) 


With  these  equations  we  can  join  together  parallel 
plate  lines  and  coaxial  lines.  We  shall  illustrate 
this  in  the  next  section. 


Fig.  1.  The  junction  of  a  parallel  plate  line  with 
a  coaxial  line.  The  structure  is  a  figure 
of  revolution  about  the  dashed  center 
line,  marked  C^. 

In  the  region  for  which  r  >  b  and  z^-h  <  z  <  z^  the 

fields  are  given  by  Eqs.  (1.4).  When  z  >  z^  and 

a  <  r  <  b  the  fields  are  given  by  Eqs.  (1.13).  The 

region  z^  -  h  <  z  <  z^,  a<r<b,  shown  dotted  in 

Fig.  1,  is  a  region  common  to  both  lines.  It  is 

bounded  by  the  two  surfaces  Z  and  2  ,  as  shown  in 

J  p  c 9 

the  figure,  by  the  plane  z  =  z^-h  and  by  the 

cylinder  r  =  a.  In  this  last  region  the  B-lines 

bend  so  that  a  B^  in  the  parallel  plate  section 

becomes  a  B  in  the  coaxial  section, 
z 


Compound  Structure 

In  this  section  we  shall  consider  a  stepped 
coaxial  line  of  n  sections  driven  by  a  parallel 

plate  line  which  in  turn  is  asymmetrically  driven  by 
a  current  source.  In  Fig.  2  we  depict  such  a 
structure  with  three  steps,  i.e.,  n  =  3.  We  shall, 
however,  keep  n  arbitrary,  as  this  causes  no 
complication.  The  coaxial  lines  and  parallel  plate 
sections  will  be  numbered  1,  2,  ...,  n  from  top  to 

bottom;  parallel  plate  line  number  n  is  the  line 

which  drives  the  whole  structure.  The  total  length 
of  the  stepped  coaxial  line  we  denote  by  L;  thus  E  = 
£,+£„+•••£.  The  quantities  a  ,  b  ,  £  (s  =  1, 

2,  •  •  • ,  n)  refer  respectively  to  the  inner  radius, 

the  outer  radius,  and  the  length  of  the  s —  coaxial 
section,  and  p  =  ^(a  +  b  ).  We  shall  denote 

s  ^  s  r  s ) 

coefficients  in  the  s —  coaxial  section  by  , 

m 
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0 


(3.4) 


and  in  the  s—  parallel  plate  section  by  A^s\ 

B^s)  . 

m 


Fig.  2.  A  stepped  coaxial  line  joined  to  a 
parallel  plate  line.  The  structure  is  a 
figure  of  revolution  about  the  dot-dashed 
center  line. 

We  specify  that  there  is  a  short  circuit  at  the 
top  of  coaxial  line  1,  i.e.,  at  the  very  top  of  the 
structure  where  z  =  L.  At  this  point  B  must 
vanish,  thus  Eq.  (1.13)  gives 


-mL/p 


E  e 
m 


1 


E  e 


mL/p, 


(3.1) 


00 

B r,  -  21/a..  +  2  2(m/a.)E  cos  m0 

u  1  l  m 

m=l 

Our  final  task  will  be  evaluate  the  quantities  E^. 

The  true  situation  we  wish  to  model  is  fed  by  a 
pair  of  rectangular  parallel  plates.  With  our 
formalism  it  is  very  difficult  to  handle  such 
shapes.  Therefore,  we  shall  make  an  assumption 
which  appears  reasonable.  Namely,  we  shall  assume 
that  we  preserve  the  essential  features  of  the 
problem  if  our  feed  line  consists  of  a  pair  of 
parallel  circular  discs,  concentric  with  the  stepped 
coaxial  line,  and  if  we  then  feed  these  discs  in  a 
realistically  asymmetric  fashion.  Accordingly  our 
feed  line  consists  of  two  discs  of  radius  b0. 

The  boundary  condition  at  the  conducting 
surfaces  of  the  parallel  plate  line  is  that 

f  4n  K  on  the  upper  plate 

B0  =  <  r  (3.5) 

(  -4n  on  the  lower  plate  , 

where  K  is  the  radial  component  of  the  surface 
current  density.  We  shall  specify  the  ultimate 
drive  of  our  system  as  a  prescribed  radial  surface 
current  density  injected  over  the  perifery  of  the 
parallel  plate  line,  r  =  b0.  Thus 

Kr  =  F(t,6)  on  r  =  bQ  ,  (3.6) 


where  E^  is  a  new  set  of  constants  (time  functions) 

introduced  for  convenience.  Let  us  now  join  this 
shorted  coaxial  line  with  the  first  parallel  plate 
section,  this  to  the  second  coaxial  section  and  so 
on.  By  induction  we  obtain  the  result 


A(n) 


=  WV 


+my 


where  F  is  a  known  function  of  t  and  6,  so  we  can 
Fourier  expand  writing 

00 

F(t,6)  =  2  F  (t)  cos  m8  .  (3.7) 

m=0  m 

According  to  Eq.  (3.5)  we  compare  Fourier 
coefficients  of  Eqs .  (3.3)  and  (3.7)  and  obtain 


B(n) 

m 


-  vvv 


Y 

n 


n 

i  «■  /p 


(3.2) 


2I/b0  =  4nFQ 


E  =  (4nb„/m)F 
m  0  m 


(3.8) 


-mV 

(bj/bjj)"1  e  " 

,2a  2mYn 
1  +  (Lj/Lq)  e 


For  our  compound  structure,  the  fields  in  the 
driving  line,  parallel  plate  section  n,  are 

»  my 

B  =  I  (m/r)E  (r/b, )m  e  n 
r  m=l  ra  1 

o  -2my 

(1  -  (b^/r)  m  e  n]  sin  mfl  ,  (3.3) 

»  my 

B„  =  2I/r  +  I  (m/r)E  (r/b.)m  e  “ 

0  m=l  ra  1 

9  -2my 

[1  +  (bj/r)2nl  e  n]  cos  m0  ,  =  0 

At  the  very  top  of  the  stepped  coaxial  structure,  at 
the  short  circuit  on  coaxial  line  1,  the  B-field 
evaluated  on  the  surface  of  the  inner  conductor  is 


If  we  substitute  this  into  Eq.  (3.4),  we  see  that 
the  total  B-field  at  the  top  of  the  stepped  coaxial 
line  is  given  by 


B  =  B0  =  2I/3l  + 


8n  £ 
al  m=l 


-my 

m  'n 


VVV  e 


2m  -2my 


1  +  (b1/bQ)  e 


cos  m8 


(3.9) 


The  infinite  sum,  which  represents  the  asymmetry  in 
the  B-field  at  the  top  of  the  compound  coaxial 
system  is  seen  to  be  small,  for  normally  Bj/b0  will 

be  quite  small  and  yn  will  be  noticeably  greater 

than  one.  Because  of  the  exponential  dependence  of 
the  coefficients  of  cos  m0  on  m,  only  the  first  few 
terms  will  be  important. 

We  conclude  this  section  with  a  few  concrete 
examples  of  asymmetric  input  current  distributions. 
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First  consider  the  case 


!K(t)  for  7t  -  0  <  6  <  n  +  9.  ,  and 

0  0 

0  otherwise 

(3.10) 


We  are  driving  the  parallel  plate  line  by  a  current 
which  is  uniform  over  a  circular  arc  of  width 
w  =  2b08.,  with  a  total  current  I  =  2b.  0.  K(t). 
Fourier  expanding  Eq.  (3.10)  gives 


F0  = 


o0k 


F  =  (-)ni  ^  sin  mO.  ,  m  t  0 
m  mn  O’ 


(3.11) 


With  these  values  for  the  F  ,  Eq.  (3.9)  becomes 

nr  ^ 


B  -  21 

Be  aa 


2  (-)""( 
n=l  \ 


1  + 


4  sin  m0_ 


mO 


0 


-my 

/i  j  •  n  id  n 
(bl/b0)  e 


,,  ..  .  2m  2mXn 

l+^j/bjj)  e 


cos  m8J 
(3.12) 


at  the  top  of  the  compound  coaxial  line. 

As  a  second  example  we  consider  a  symmetric  two 
sided  feed.  Thus  we  set, 


F(t,0)  = 


yc(t)  for  n  -  0q  <  0  <  71  +0q  and  for 


90  <  0  <  -eo  >  and 

0  otherwise 


(3.13) 


Now  our  Fourier  coefficients  are 


F 

0  7t 


F  =  [1  +  (-)"1]  sin  ra0  ,  m  *  0 

m  m7i  0 


(3.14) 


For  even  values  of  m  this  is  identical  with 
Eqs.  (3.11)  whereas  for  odd  values  of  m,  F 
vanishes.  Thus,  for  the  field  at  the  top  of  the 

coaxial  line  we  recover  Eq.  (3.12)  modified  solely 
by  deleting  from  the  summation  all  terms  for  which  m 
is  odd. 

Generally,  for  q  equally  spaced  feed  strips  we 
simply  delete  from  Eq.  (3.12)  those  terms  for  which 
m  is  not  divisible  by  q. 


Discussion  of  the  Asymmetry 

As  an  illustration  of  the  foregoing  analysis, 
let  us  consider  a  numerical  example,  which  is  an 
early  design  for  one  of  our  foil  implosion 
experiments.  The  foil  is  the  very  top  section 
(line  1)  of  Fig.  1.  The  field  is  given  by 
Eq.  (3.12);  B  and  B  vanish  here.  Our  stepped 
coaxial  line  consists  of  three  sections ,  whose 
dimensions  are  given  in  Table  1. 

Table  1. 

Dimensions,  in  cm.,  of  the  Stepped  Coaxial  Line 


a , 

1 

b. 

1 

pi 

1 

£./p. 

1  *1 

4.2 

4.7 

4.45 

5.5 

1 . 2340 

8.2 

8.7 

8.45 

3.4 

0.4024 

11.0 

11.5 

11.25 

4.8 

0.4267 

For  this  line  n  =  3  and  h  =  0.5  cm.  From  the  table 

-y 

we  get  Ya  s  ~  2.063  and  e  n  =  0.1271 

The  two  feed  plates  are  square,  76.2  cm 
(30  in.)  on  a  side  and  are  symmetrically  fed  over 
the  full  width  of  two  opposite  sides.  Thus  0^,  the 

half  angle  of  one  of  the  feed  strips,  is  rt/4,  whence 
sin  0.  sin  20- 

— x — -  =  0.9003  and  =  0.6366  .  (4.1) 

60  290 


We  must,  for  our  model,  replace  the  square  feed 

plates  by  circular  discs.  We  conservatively  take 
bg  =  38.1  cm.,  i.e.,  we  replace  the  squares  by 

inscribed  circles.  This  gives  b^/b^  =  0.1234  whence 

X  =  (b./b  )e  n  =  1.57  x  io  , 

2-4  <4‘2> 

X  =  2.46  x  io  . 


For  such  a  geometry  we  see  that  all  the  denomenators 
in  Eq.  (3.12)  can  be  replaced  by  unity  and  the  first 
three  terms  become 


21  4  sin  60 

B0  =  ^  f1  *  <— 0^X  cos  6 


4  sin  20 

+  ( - 20 - ^)X  cos  26  •••  J 


(4.3) 


if  we  drive  the  parallel  plates  from  one  side  only. 

The  importance  of  the  field  asymmetry  is  that 
the  consequent  asymmetry  of  the  magnetic  pressure 
will  drive  a  Rayleigh-Taylor  instability.  The 

2 

driving  pressure  is  P  =  Bg/87t.  Using  the  above 
numerical  values  we  would  have 

I2 

P  =  — — -r  {1  -  0.1131  cos  0 

27ta1Z 


+  (2.852  x  io'3)  cos  26  •••} 

(4.4) 

if  the  line  were  driven  from  one  side  only. 

For  the  balanced,  two  sided  drive,  we  must  omit 
the  terms  of  Eq.  (4.3)  in  cos  6,  cos  30  etc.  In 
this  case  we  obtain 


P  =  — —“j  {1  +  (1.253  x  10*3)  cos  20  +  •••} 

2nal  (4.5) 

From  the  value  of  X  in  our  case,  Eq.  (4.2),  it  is 
clear  that  the  further  terms  ”in  Eqs.  (4.4)  and 

3  4 

(4.5),  proportional  to  X  cos  30,  X  cos  40  etc., 
are  completely  negligible. 

Having  calculated  the  pressure  asymmetry,  Eq. 
(4.4)  or  (4.5),  we  have  the  input  data  for  our 
principal  problem,  namely  describing  the  resulting 
Rayleigh-Taylor  instability.  This,  however,  is  the 
subject  for  a  different  paper. 
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